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We discuss laser dressed dipolar and Van der Waals interactions between atoms and polar 
molecules, so that a cold atomic gas with laser admixed Rydberg levels acts as a designed reservoir 
for both elastic and inelastic collisional processes. The elastic scattering channel is characterized by 
large elastic scattering cross sections and repulsive shields to protect from close encounter collisions. 
In addition, we discuss a dtssipative (inelastic) collision where a spontaneously emitted photon car- 
ries away (kinetic) energy of the collision partners, thus providing a significant energy loss in a single 
collision. This leads to the scenario of rapid thermalization and cooling of a molecule in the mK 
down to the /jK regime by cold atoms. 

PACS numbers: 34.20. Gj, 37.10. Mn, 32.80. Ee 



There is at present significant interest in preparing 
and manipulating cold samples of molecules [THS]. A 
promising avenue towards this goal seems to employ the 
ubiquitous ultracold atomic gases as cold reservoirs, and 
to study mixtures of atomic and molecular gases, where 
molecules and atoms interact via collisional processes [4] . 
Given the well developed tools in manipulating atoms 
with external electromagnetic fields [S], it is natural to 
ask whether we can "design" these atom-molecule inter- 
actions, thus effectively engineering an atomic reservoir 
with desired (collisional) properties. Below we will de- 
scribe a specific scenario of engineered elastic and inelas- 
tic collisions involving laser-dressed atoms and ground 
state molecules with remarkable, and potentially useful 
properties. This includes (i) strong repulsive shields to 
protect from inelastic collisions and chemical reactions 
and (ii) exceedingly large scattering cross sections for 
elastic scattering between the atom and the molecule. 
The relevant energy (temperature) range includes sev- 
eral mK down to fiK. Equally important, we will show 
that (iii) we can design a "dissipative collision" where a 
spontaneously emitted photon carries away (kinetic) en- 
ergy of the collision partners, thus providing a significant 
energy loss in a single collision what could be called "col- 
lisional Sisyphus" effect [5], in analogy to Sisyphus laser 
cooling of single atoms in external trapping potentials 
[7]. This suggests rapid thermalization and cooling of a 
molecule by the cold atom reservoir. 

The atomic and molecular level scheme, and the colli- 
sional process we have in mind are illustrated in Figs.[l|a) 
and (b), respectively. The basic ingredient is the long 
range dipolar interaction between molecules in the rovi- 
brational ground state and laser excited Rydberg atoms. 
The Rydberg state |r) is chosen so that its electric 
dipole transitions to neighboring states approximately 



[*] These authors contributed equally to this work. 




FIG. 1: (a) Energy levels of a laser excited atom and a 
rotational spectrum of a polar molecule. The Rydberg state 
|r) interacts with the molecule via a dipole-dipole interaction 
Vdd (see text), (b) Born-Oppenheimer (BO) potentials for 
the laser dressed atom + molecule complex. We consider 
a dissipative collision, where (1) the particles collide on the 
potential curve V±(r) with the atom in \gi), climb the "blue 
shield" step at r c , and (2) are quenched to the potential V<i(r) 
with atom in <?2 • The dominant atomic state is indicated with 
the molecule in its ground state, (c) Decay rate 71 (r) of the 
BO potential (see text). 



matches the rotational excitation spectrum of the polar 
molecule with frequencies in the microwave regime (cf. 
Fig. [T^i), implying a near resonant exchange of molec- 
ular and atomic excitations. We focus below on the 
conceptually simplest configuration, where this interac- 
tion reduces to a large repulsive and isotropic Van der 
Waals interaction, V v dw( r ) — C-e/r 6 , a situation analo- 
gous to the large Rydberg-Rydberg interactions under- 
lying the dipole blockade mechanism and the formation 
of superatoms [SHTTj. This interaction is admixed to 
the atomic ground state \gi) with a blue detuned laser 
A r > 0, thus providing an effective interaction between 
the ground state atoms and molecules. The relevant 
Born-Oppenheimer (BO) potential for the laser dressed 
complex is sketched in Fig.[T|b) as V\{r) (see below), and 
defines a collision channel for a molecule and atom ini- 
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tially in \gi). In this collision the particles moving adia- 
batically on Vi (r) will encounter a steep "blue shield" po- 
tential !l2j at a distance Kdw( r c) = A r , where typically 
r c > 100 nm. By adding a second laser which couples the 
long-lived Rydberg state |r) down to a low- lying short 
lived excited state |e) with a detuning A e = S + A r > 
(see Fig. [iji), we can add a plateau for r < r' c to the 
adiabatic potential V\ (r), so that atoms are efficiently 
quenched to the ground state \g 2 ) with potential V 2 (r) 
according to the rate 7i(r) in Fig. [ljc). 

This leads to the following overall picture of collisions 
illustrated in Fig. [ljb): (i) for a kinetic energy of rel- 
ative motion less than the potential step in V\(r), i.e. 
-Eidn < <5, we have an elastic collision from an effective 
hard core potential with (large) radius r c ; (ii) in a colli- 
sion with .Ekin > $ the particles will climb the potential 
step entering the flat dissipative region, which acts as a 
"trap door" so that in a single collision the kinetic energy 
~ 6 is carried away by the spontaneous photon with the 
atom being left in 152)- Below we will work out a quanti- 
tative description of these collisional processes, and argue 
that they can occur with high fidelity. An essential argu- 
ment is that during the collision the particles never enter 
the small distance regime [shaded region in Fig. |TJb)] , 
and thus the collisional dynamics does not couple signif- 
icantly to other channels. The above collision cycle can 
be repeated by pumping atoms from \g 2 ) back to \g±) so 
that a significant amount of energy can be lost in a few 
collisions. Besides, due to the large collision cross sec- 
tions for elastic processes there is efficient thermalization 
of the molecules and atoms (sympathetic cooling). 

Master equation: The dynamics of the atom and 
molecule for a dissipative collision is described by a mas- 
ter equation, p = —i [H, p] + Cp, with the RHS as sum 
of a Hamiltonian and dissipative part, and p the reduced 
system density operator after tracing over the vacuum 
modes of the radiation field. Such an equation is readily 
written down as an extension of the familiar master equa- 
tions of laser cooling for atoms by including the molec- 
ular dynamics and the atom-molecule interactions. We 
neglect, however, recoil kicks from laser absorption and 
spontaneous emission, and Doppler shifts, as they pro- 
vide only small corrections to our collisional dynamics. 

The Hamiltonian has the form H — T + Hj (r) , where 
f = P 2 /2A/ + p 2 /2/i is the kinetic energy with P (p) the 
center of mass (relative) momentum and M (p) the total 
(reduced) mass, and Hj(r) = H m + Hqa + Vdd(r) is the 
Hamiltonian for the internal degrees of freedom as sum 
of a molecular and atomic Hamiltonian, and the dipole- 
dipole interaction. For the molecule we assume a rigid 
rotor Hamiltonian i?oM = -BN 2 with B the rotational 
constant, and N the angular momentum. For the atomic 
Hamiltonian we write in the rotating wave approximation 
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with notation &ij = for the atomic transition oper- 
ators, and atomic states according to Fig. 1(a). Here we 
consider the conceptually simplest situation where |r) = 
1 71, s) and \r') = \n — l,p), with n the principal quan- 
tum number and s and p the orbital angular momentum 
quantum number [T3]. By A r (5) and Q r (fi e ) we denote 
the detuning and Rabi frequency of the exciting (quench- 
ing) lasers, respectively, and by —A r , the energy of the 
\r') state. For an energy mismatch E = E a — E m <C 
mm{E a , E m } between the atomic Rydberg states, E a = 
E r — E rh and the rotational splitting of the molecule, 
E m = 2J5, the interaction is dominated by the dipole- 
dipole interactions between the two channels |r)|0) and 
|r')|l). The interaction Hamiltonian between a Rydberg 
atom and the molecule separated by a distance r = rr is 
V dd = [d. r -d m -3(d r -f)(d m -r)]/r 3 (|r)(r'|®|0)(l|+h.c.)for 
distances larger than the size of the Rydberg atom, r > 
r s ~ n 2 ao, where ao is Bohr's radius. Here, d r = (r|d|r') 
and d TO = (0|d|l) are the atomic and molecular tran- 
sition dipole moments, respectively, with d the dipole 
operator. For distances r > r CQ = (d r d m / ' E$) 1 ^ , we 
can adiabatically eliminate and obtain the effec- 

tive interaction between \r) and |0), which is a repul- 
sive and isotropic van-der-Waals interaction Ce/r 6 in 
three- dimensions with strength Cg = 2<i 2 e? 2 ra /3£t)- Fi- 
nally, the Liouvillian C in the master equation describes 
dissipative processes due to spontaneous emission. We 
write L = C e + C p + Cb, where C e = 7 e 2?[<T eg2 ] and 
C p — j%D [o^g 2gi ] account for spontaneous decay from |e) 
to I g 2 ) and re-pumping from ^2) to |<7i), respectively, 
with Lindblad term T>[a]p = apa^ — a^ap/2 — pa^ a/2. 
The third term Cbp describes undesired decays, includ- 
ing in particular spontaneous emission from the Rydberg 
state, as discussed below. 

Born-Oppenheimer approximation: We proceed by 
identifying the BO potentials of the dressed atom- 
molecule complex as eigenvalues of the internal Hamilto- 
nian Hj(r)\i(r)) — Vi(r)\i(r)) depending parametrically 
on r [compare Fig. 0b)], which, in an adiabatic approxi- 
mation, provide effective interaction potentials for atoms 
and molecules. In particular, the dressed groundstate po- 
tential V\(r) corresponds to the BO energy surface that 
asymptotically connects to the ground state of the atom 
at large distances, i.e., |l(r — > 00)) ~ |ffi,0). There, 
atom and molecule are essentially non-interacting. As ex- 
plained above, the step-like character is obtained in com- 
bination with laser dressing on two internal atomic tran- 
sitions: (i) by coupling \gi) with |r) in the weak-dressing 
regime Q r /A r < 1 and for blue detuning A r > 0, V\(r) 
becomes approximately V\ (r) ~ Ca/r 6 for distances r < 
r c , with r ~ r c = (Cs/Ar) 1 ^ a resonant Condon point 
with typical values in the hundreds of nm. This design of 
interactions is similar to blue-shielding techniques with 
cold atoms and molecules, however it exploits repulsive 
vdW-interactions and thus works in three-dimensions. 
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The dominant contribution to the 1 1 (r)) is now \r, 0), 
with r r ~ n 3 the lifetime of |r), e.g., in the hundreds of 
/zs regime for n ~ 80 [ID] , (ii) A second Condon point 
can be engineered at distances r£ = (Cg/A,,) 1 / 6 < r c by 
weakly admixing \r) with the low-energy excited state |e), 
using laser light with fl e /A e < 1 and A e > A r . Here we 
assume that |e) interacts only weakly with the molecule 
and thus V\{r) becomes essentially flat for r < r' c . Popu- 
lation in |e) quickly decays to a second groundstate \g 2 ) at 
a rate 7 e ~ MHz. This makes the decay rate from 1 1 (r)) 
strongly position-dependent as 71 (r) = 7 e |(l(r)|e,0)| 2 , 
see sketch in Fig. [ljc). The BO potential V2(r) with 
|2(r)) ~ |g 2 ,0) is essentially flat, Fig.gb). 

Different BO potentials are coupled via residual non- 
adiabatic transitions at r c and r' c . In particular, popula- 
tion transfer at r c from |l(r)) to the BO eigenstate that 
connects to |r, 0) for r ^> r c could induce significant heat- 
ing and losses. An estimate of the non-adiabatic transi- 
tion probability can be computed for b -C r c within a ID 
Landau Zener model as Plz = cxp(— 27ril 2 ,/(au)), with 
a is the difference of the gradient of the bare potentials 
at r c . This shows that for any given velocity v (in rela- 
tive coordinates) non-adiabatic transitions can be always 
suppressed by increasing fi r . Full 3D computations of the 
non-adiabatic transition probabilities in the semiclassical 
limit confirm these predictions, see Appendix. Since |e) 
decays to \gz), diabatic transitions at r' c from 1 1 (r)) to 
the BO eigenstate which adiabatically connects to |e, 0) 
for r 3> r' c are allowed in our scheme. 

Additional diabatic crossings with potential surfaces 
involving different Rydberg states as well as attractive 
resonant dipole-dipole interactions lead to collisional two- 
body losses for distances r < r Co . Moreover, interactions 
between the Rydberg-electron and the molecule play a 
significant role for r < r s [14!. These effects are sup- 
pressed by a "blue-shield" at r"> max{r s ,r CQ } confining 
particles' motion to r > r" [To] . 

Reservoir engineering and molecular cooling: In our 
scheme, we consider hot molecules undergoing a few scat- 
tering processes with cold, interacting, Rydberg-dressed 
atoms, with lifetime ~ r r (f2 r /2A r )~ 2 . Cooling of the 
molecules comes as a combination of sympathetic cooling 
with atoms with large elastic cross sections a ~ 7rr 2 as 
well as photon-assisted controlled inelastic interactions in 
a timescale r c < to avoid spontaneous emission from 
\r) and collisional losses with Rydberg-excited atoms [ID] . 

The basic scheme of photon-assisted inelastic collisions 
can be understood for just an atom and a molecule ini- 
tially interacting via the BO potential V\ (r) . It comprises 
two steps: Firstly, spontaneous emission from |e) couples 
I l(r)) and |2(r)), according to the spatially-dependent 
rate 7i(?"), removing an amount of energy < 5] secondly, 
a weak re-pumping laser can transfer population back 
from |2(r)) to |l(r)), thus closing the cooling cycle. By fo- 
cussing on V\{r) and V^r) only and neglecting for a mo- 
ment unwanted effects described by Cbp, within the secu- 
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FIG. 2: Dissipative collisions: (a) Energy loss Ei oss per colli- 
sion vs impact parameter b and initial relative kinetic energy 
Skin- (b) Distribution / of final relative kinetic energies ££ in 
after a single collision for different b, for an initial Boltzmann 
distribution with ksT/S ~ 0.5. (c) Average final kinetic en- 
ergy (E' kin ) of the corresponding distribution after a single 
collision, (d) E k i n (solid lines, left axis) and the average rel- 
ative kinetic energy (Skin) (dotted line, right axis) vs time 
with r = povy, ss 2-7T x 2.3 kHz the collision rate (see text). 
Parameters: d m = 7 Debye, d r = 4400 Debye, Eq = 2n x 2.5 
GHz, A r = 2tt x 60 MHz, fl r /A r = 0.42, S/A r = 0.33, 
QJA r = 0.25, r c « 218 nm, r'Jr c w 0.95 , and r"/r c « 0.81. 



lar approximation the dissipative collisional dynamics in 
the relative-coordinate frame can be described semiclas- 
sically by two coupled Liouville-equations 
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Here /i(r, p, t) is the Liouville density accounting for the 
phase-space distribution of the atom-molecule system in 
state \i(r)) £ {l,2},i 7^ j). The first term in the 
RHS is the interaction force, proportional to the gradi- 
ent of the BO potentials discussed above. The second and 
third terms are the spatially dependent decay rate 71 (r) 
and the re-pumping rate 72, coupling the two equations. 
The step-like shape of 71 (r) sketched in Fig. [ljc) reflects 
the one of Vi(r), such that 71 (r) ~ j e for r" < r < r' c and 
71 (r) ~ otherwise. As a result, for incoming relative 
kinetic energies E k [ n < 5, particles are reflected elasti- 
cally at r rs r c , while particles with E'kin > $ can reach 
the region r < r' c with a velocity v' = ^/2(i?kin — S)/fi 
and undergo photon-assisted inelastic collisions. For any 
given v', spontaneous emission from |e) to \gz) [and thus 
population transfer from | l(r)) to |2(r))] can be made to 
occur deterministically in a region of length d = r' c — r" , 
by choosing d such that jed/v' > 1. This removes an 
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energy of order E\ oss ~ S in every single collision, with 
8 as large as mK, as shown below. We ensure that the 
inverse re-pumping from |2(r)) to 1 1 (r)) takes place for 
distances r > r c by requiring 72<i/V C 1. In addition, 
we choose 7 e [f2 r f2 e /(4A r i5)] 2 < 72 to ensure that the ef- 
fective Raman transfer rate of population from \g\) to 
|g 2 ) via \r) is small, and thus the atomic population is in 
I <7i) at r 3> r c . These requirements can be satisfied for 
realistic atom/molecule configurations, as shown below. 

We investigate numerically this dissipative scheme by 
performing molecular dynamics simulations of the colli- 
sion of an atom and a molecule, based on Eq. ([!]). The 
extension to the case of several atoms and molecules is 
straightforward. In a semiclassical approximation, the 
mean energy loss -Ei OS s in a collision in 3D is computed 
as E loss = f[Vi(r) - V 2 (r)]7i(r c i)p(r cl )cit, where r d de- 
notes the classical trajectory of the atom/molecule colli- 
sion, and p(r c i) is the probability that the atom decays 
at a given position r c i, with p(t) — — 7i(r c i)p(£). As an 
example, in the calculations we consider a NaH molecule 
(d m ~ 7 Debye) and a Cs atom with Rydberg states 
|r) = |46s) and \r') = |45,p), respectively, withd r ps 4400 
Debye and Eq = 2ir x 2.5 GHz. The laser parameters 
are A r = 2tt x 60 MHz, Vt r /A r = 0.42, 8/A r = 0.33, 
fl e /A r = 0.25, so that r c « 218 nm, r' c « 208 nm, and 
r" w 177 nm [T7]. Figure^ a) shows the computed Ei oss 
as a function of the initial relative kinetic energy Ekin 
and the impact parameter b. For b > r c the collision is 
essentially elastic, as expected. However, for b < r c and 
energies E^ n > 5 the molecule is able to climb the poten- 
tial step (5 of V\ at r ~ r c , thus undergoing deterministic 
decay to |2(r)). 

The effects of this dissipative collisional cooling on a 
molecular gas with an initial thermal distribution with 
average temperature T = 0.5 mK where 8 ~ 2k B T with 
ks Boltzmann's constant is shown in panel (b). For 
each fixed value of b we perform N ~ 10 5 computa- 
tions of the collision dynamics, by randomly generating 
a sample of initial kinetic energies -Ekinj according to a 
Boltzmann distribution. The figure shows the final pop- 
ulation distribution / as a function of the final energy 
E' kin = Skin — ^loss, for fixed values of b, with laser pa- 
rameters as in Fig. [5Ja). We find that for impact pa- 
rameters b > r c the distribution / in relative coordinates 
is largely unaffected by the collision (case b — r c in the 
figure). For b < r c , however, all population with initial 
energy Ekin > 8 is shifted by an amount ~ 8 towards 
lower energies. The corresponding average final E^ in is 
shown in Fig. [2jc) as a function of b. For heads-on colli- 
sions with 6 = approximately 50 % of the initial kinetic 
energy is removed after a single collision. 

For given 8, E^ n and atoms at rest, the energy loss 
rate is estimated as —dE^/dt = p{2E^\ n / ' p) 1 / 2 J r (E\^ n ), 
with -F(iSkin) = Jq" E\ oss (b, Eki n )2nbdb and p the atomic 
gas density. The latter is limited by atom-atom inter- 
actions of the form V aa ~ (il r /(2A r )) 4: V rr for atomic 



distances r > Pmax 3 , with V rr = Ce/r 6 the vdW interac- 
tion between Rydberg states and pmlx = [Ce/(2A r )} 1 ^ 6 
a resonant Condon radius (p m ax 3 — 1.7 /im for the pa- 
rameters above). Figure J2jd) shows the presence of two 
cooling timescales (solid lines): For £kin > 8, cooling 
of an energy ~ 8 is achieved on a fast timescale of a 
few Ft, with T — pav^ w 2tt x 2.3 kHz the collision 
rate, for p = (2 ^m) -3 , r c w 218 nm, T = 0.5 mK, 
Wju = ^ik B T/ p ss 0.78 m/s. For i? kin < <5 cooling 
proceeds slowly, in accordance with the small 71 (r), for 
r > r c . The same qualitative behavior is found in the av- 
erage kinetic energy (dots). The lifetime of the dressed 
state is here s=s 1 ms (r r ss 45 /xs), and thus for the pa- 
rameters above more than 10 collisions are allowed while 
cooling. We note that 8 can be dynamically reduced in 
experiments. 

In the lab frame, a molecule loses its energy due to 
a combination of both collisional dissipative cooling and 
sympathetic cooling. The dominant effect depends on 
the mass ratio m^/mjvf- For an atom initially at rest, an 
analytic estimate for the atomic and molecular velocities 
v' A and v' M after the collision can be obtained from a 
simplified model where the total energy is reduced by 
8 whenever p,v 2 M /2 > 8, with vm the initial molecular 
velocity, and is conserved otherwise, as 

v' M = V(l - m A /m M \]l - 28/(nv 2 M )) 

v' A = V(l + y/l-28/(nv 2 M )). (2) 

Here, V = ttimVm/M is the center of mass velocity. 
Figure [3] shows the result of molecular dynamics simu- 
lations where we study the energy loss of the molecule 

^loss = ^kir? ~ 0-/2) m M'v'M f° r different mass ratios 
and laser parameters as in Fig. [2] In the figure, the 
dashed and continuous lines correspond to pure sympa- 
thetic cooling and the predictions of Eqs. pj), respec- 
tively, while squares and dots are numerical results for 
different values of , averaged over 200 simulations. 
For rriA < tum the dominant energy loss mechanism is 
sympathetic cooling. However, for m A > mjj the en- 
ergy loss of the molecule is mainly caused by dissipative 
collisional cooling, and is of the order of 8, as expected. 
The effective atomic mass may be tuned using external 
confining potentials, e.g., optical lattices. For example, 
atoms can be confined in an optical trap with depth 
Vtr 2m aV 2 , which for NaH with vm = \/3ksT '/tum 
and T = 0.5 mK and Cs atoms with rriA/mM ~ 5.5 
implies V tt > 0.78k B T w 0.4 mK. 

In conclusion, we have discussed a scenario where 
a molecule scatters successively from cold (stationary) 
atoms in designed elastic and inelastic processes. In this 
situation reminiscent of a "microscopic version of a pin- 
ball machine", inelastic scattering events are associated 
with the emission of a photon implying a "collisional Sisy- 
phus" cooling. While we focused on the simplest possi- 
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FIG. 3: Energy loss vs mass ratio mAftnu- The dashed 
red and solid black lines are analytic results for a model of 
sympathetic cooling only (S = 0) and finite 5, respectively, 
see Eqs. |2| and text. Blue dots and squares are averages 
over 200 runs of molecular dynamics simulations for finite S, 
for laser parameters as in Fig. [2] 

ble setup based on Van der Waals interactions, variants 
based on, e.g. dipole-dipole interactions and low dimen- 
sional trapping geometries seem possible. We will inves- 
tigate the role of many-atom interactions in the dynamics 
of the gas in future work. 

Note added: In the final stages of work we became 
aware of S.D. Huber and H.P. Biichler's proposal for 
Doppler cooling of polar molecules, where atomic Ryd- 
berg excitations serve as a bath for rotational molecular 
excitations [15] , 

We thank F. Herrera for discussions, and S.D. Huber 
and H.P. Buchler for sharing with us their work before 
publication. Work supported by the Austrian Science 
Fund, EU grants AQUTE, COHERENCE and NAME- 
QUAM, and by MURI, AFOSR, and EOARD. 

Appendix 

Non-adiabatic transition: We calculate the classical 
trajectory by fir'd = — VVi(r c i). Plugging the trajectory 
into the Schrodinger equations governing the dynamics of 
internal states and calculating the transition probabilities 
after the molecule has reached the flat region. The results 
are shown in Fig. [4j with laser parameters as in Fig. [2] p\ 
is the computed transition probability from 1 1 (r) ) to the 
BO eigenstate that connect to |r, 0) for r ^> r c , which is 
on the order of 10 -2 for large kinetic energies -Ekin < 3<5. 
P2 is the transition probability from | l(r)) to the BO 
eigenstate that connect to |e, 0) for r ^> r' c , which is 
tolerant as discussed. P3 is the non-adiabatic transition 
probability in |2(r)) at r", which is calculated in a simi- 
lar way. Note that non-adiabatic transition probability in 
|l(r)) at r" is not important, since spontaneous decay al- 
most takes place deterministically. All the non-adiabatic 
transitions can be further suppressed by increasing the 
Rabi frequency. 



(a) (b) 




FIG. 4: Diabatic transitions between different BO eigenstates. 

(a) pi versus impact parameter b and initial kinetic energy. 

(b) Transition probabilities for 6 = 0. 
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